Toeplitz matrices: spectral properties and preconditioning in the CG method
نویسنده
چکیده
We consider multilevel Toeplitz matrices Tn(f) generated by Lebesgue integrable functionsf defined over I, I = [−π, π), d ≥ 1. We are interested in the solution of linear systems withcoefficient matrix Tn(f) when the size of Tn(f) is large. Therefore the use of iterative methodsis recommended for computational and numerical stability reasons. In this note we focus ourattention on the (preconditioned) conjugate gradient (P)CG method and on the case where thesymbol f is known and univariate (d = 1): the second section treat spectral properties of Toeplitzmatrices Tn(f); the third deals with the spectral behavior of T−1n (g)Tn(f) and the fourth with theband Toeplitz preconditioning; in the fifth section we consider the matrix algebra preconditioningthrough the Korovkin theory. Then in the sixth section we study the multilevel case d > 1 byemphasizing the results that have a plain generalization (those in the Sections 2, 3, and 4) andthe results which strongly depend on the number d of levels (those in Section 5): in particularthe quality of the matrix algebra preconditioners (circulants, trigonometric algebras, Hartley etc.)deteriorates sensibly as d increases.A section of conclusive remarks and two appendices treating the theory of the (P)CG methodand spectral distributional results of structured matrix sequences.key words Linear system, conjugate gradient method, Toeplitz matrix, structured matrix, precon-ditioner.
منابع مشابه
On the Spectrum of a Family of Preconditioned Block Toeplitz Matrices
Abstract. Research on preconditioning Toeplitz matrices with circulant matrices has been active recently. The preconditioning technique can be easily generalized to block Toeplitz matrices. That is, for a block Toeplitz matrix T consisting of N N blocks with M M elements per block, a block circulant matrix R is used with the same block structure as its preconditioner. In this research, the spec...
متن کاملA Preconditioned MINRES Method for Nonsymmetric Toeplitz Matrices
Circulant preconditioners that cluster eigenvalues are well established for linear systems involving symmetric positive definite Toeplitz matrices. For these preconditioners rapid convergence of the preconditioned conjugate gradient method is guaranteed. Since circulant preconditioners can be applied quickly using the fast Fourier transform, preconditioned CG with circulant preconditioning is e...
متن کاملPreconditioning Methods for Shift-Variant Image Reconstruction
Preconditioning methods can accelerate the convergence of gradient-based iterative methods for tomographic image reconstruction and image restoration. Circulant preconditioners have been used extensively for shiftinvariant problems. Diagonal preconditioners offer some improvement in convergence rate, but do not incorporate the structure of the Hessian matrices in imaging problems. For inverse p...
متن کاملSpectral analysis and preconditioning techniques for radial basis function collocation matrices
Meshless collocation methods based on radial basis functions lead to structured linear systems which, for equispaced grid points, have almost a multilevel Toeplitz structure. In particular, if we consider partial differential equations (PDEs) in two dimensions then we find almost (up to a “low-rank” correction given by the boundary conditions) two-level Toeplitz matrices, i.e., block Toeplitz w...
متن کاملSpectral behavior of preconditioned non-Hermitian multilevel block Toeplitz matrices with matrix-valued symbol
This note is devoted to preconditioning strategies for non-Hermitian multilevel block Toeplitz linear systems associated with a multivariate Lebesgue integrable matrix-valued symbol. In particular, we consider special preconditioned matrices, where the preconditioner has a band multilevel block Toeplitz structure, and we complement known results on the localization of the spectrum with global d...
متن کامل